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Abstract 

The standard interference functions introduced by Yates have been very influential on the analysis 
and design of distributed power control laws. While powerful and versatile, the framework has some 
drawbacks: the existence of fixed-points has to be established separately, and no guarantees are given on 
the rate of convergence of the iterates. This paper introduces contractive interference functions, a slight 
reformulation of the standard interference functions that guarantees the existence and uniqueness of 
fixed-points along with linear convergence of iterates. We show that many power control laws from the 
literature are contractive and derive, sometimes for the first time, analytical convergence rate estimates 
for these algorithms. We also prove that contractive interference functions converge when executed 
totally asynchronously and, under the assumption that the communication delay is bounded, derive an 
explicit bound on the convergence time penalty due to increased delay. Finally, we demonstrate that 
although standard interference functions are, in general, not contractive, they are all para-contractions 
with respect to a certain metric. Similar results for two-sided scalable interference functions are also 
derived. 



I. Introduction 

Distributed power control (DPC) algorithms such as and Q have had an enormous 
influence on modern wireless systems. The basic algorithm for adjusting transmit powers to 
meet predefined Signal-to-Interference-and-Noise-Ratio (SINR) targets can be written as a linear 
iteration and has been thoroughly analyzed using tools from linear algebra and positive linear 
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systems. In particular, when the SINR targets are feasible, the algorithm converges to a unique 
fixed-point at a linear rate, i.e., the distance between the iterates and the optimal power allocation 
decays exponentially. These results can be derived using Perron-Frobenius theory for positive 
matrices or, alternatively, by showing that the linear iteration is a contraction mapping in a 
weighted maximum norm (see for example [3J). 

An elegant axiomatic framework for studying more general power control iterations was 
proposed by Yates H). The so-called standard interference functions include the linear iterations, 
and several important nonlinear power control laws. Various extensions of the basic framework 
have been proposed in the literature with the most prominent those by Sung and Leung Q 
and Schubert and Boche [|6j|. While several results exist for synchronous and asynchronous 
convergence of interference function iterations, very few results on the convergence rates of 
such algorithms have appeared in the literature (see for example |I3-|l9l for exceptions). The 
current proofs are tailor-made and the link to contraction mappings, that has been so powerful 
in the analysis of the linear iterations, is disturbingly absent. This paper tries to fill this gap. 

Contrary to claims in the literature, we demonstrate that interference functions are, in general, 
not necessarily contraction mappings. However, we show that a slight modification of the 
scalability axiom of standard interference function allows to guarantee contractivity of the 
iterations and hence unique fixed-points and linear convergence rates. This condition is satisfied 
by the basic DPC algorithm and allows to recover the same convergence rate that comes out of 
a tailored analysis. It also allows to estimate the convergence rate of the other power control 
schemes considered by Yates, as well as the utility-based power control scheme developed in [lOJ. 
We also demonstrate how a logarithmic change-of-variables render interference functions para- 
contractions. Interestingly, this is the same change of variables that has been used extensively 
in resource allocation for interference-limited systems {e.g., ifTTIl . Furthermore, we show 

that our conditions are also satisfied for the two-sided scalable interference functions introduced 
in [m, and introduce conditions that guarantee that two-sided scalable interference functions 
define contraction mappings, and hence have unique fixed-points and linear convergence rates. 
Finally, we discuss how asynchronous convergence can be established in our framework. 
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A. Notation 

Throughout the paper, vectors are written in bold lower case letters and matrices in upper case 
letters. The set of non-negative real numbers is denoted by IR+. Xi denotes the i*^ component 
of a vector x, and the notation x > means that all of the components of x are greater than or 
equal to zero. The inequality x > y implies that Xi > Ui for all components i. We use e'' and 
ln(x) to denote component-wise exponential and logarithm of the entries x, respectively. The 
matrix A is said to be Hurwitz if all its eigenvalues have negative real parts. It is Metzler if all 
its off-diagonal entries are non-negative. Given a vector v > 0, || ■ ||^ stands for the weighted 
maximum norm, i.e., ||x||JL, = max 



denoted by || • ||^, defined by 



The vector norm || ■ ||^ induces a matrix norm, also 



~ x^O X v 



When Vi = 1 for all i, we suppress the superscript v. The spectral radius of a matrix A is the 

largest magnitude of the eigenvalues of A and is denoted by p{A). A sequence {x„} G is 

said to converge linearly to x* if there exists a constant c G (0,1) such that 

||a;(n+ 1) - x^ll 
lim 11 , . = c, 

n^oo \\X[n) — X \\ 

where || ■ || is some norm on M^^. If c = 1, then the sequence is said to converge sub-linearly. 

II. Fixed-point theory and interference functions 
A. Fixed-points, contractions and para-contractions 
We consider iterative algorithms on the form 

x(n + l) = T(x(n)), n = 0,l,2,..., (1) 

where T is a mapping from a subset X of into itself. A vector x* is called a fixed point of 
T if T(x*) = X*. If T is continuous at x* and the sequence {x(n)} converges to x*, then x* is 
a fixed point of T [[T3l Chapter 3]. Therefore, the iteration ^ can be viewed as an algorithm for 
finding such a fixed point. T is called a contraction mapping, if it has the following property 

||T(x)-T(y)|| <c||x-y|| ,Wx,y e X, 

where || ■ || is some norm on X, and c G [0, 1). The following proposition shows that contraction 
mappings have unique fixed points and linear convergence rates. 
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Proposition 1 (Convergence of Contracting Iterations |[13l Chapter 3]) If T : X ^ X is a 

contraction mapping and that X is a closed subset of M^, then: 

• (Existence and Uniqueness of Fixed Points) The mapping T has a unique fixed point x* G X. 

• (Linear Convergence) For every initial vector x(0) G X, the sequence {x(?7,)} generated 
by x(n + 1) = T(x(n)) converges to x* linearly. In particular, 

||x(n) -x'^ll < c"||x(0) -x^ll . 

An operator T on X is called para-contraction if 

||T(x)-T(y)|| < ||x-y|| , for all x 7^ y . 

Para-contractions have at most one fixed point and, in contrast to contractions, may not have a 
fixed point. As an example, consider the para-contracting function T(x) = x + e~^' in [0,oo). 
It is easily seen that T has no fixed point. The following theorem summarizes properties of 
para-contractions . 

Proposition 2 ( |T4l) If T : X ^ X is a para-contraction, then: 

• IfT has a fixed point x*, then that fixed point is unique; moreover 

• If X is a finite -dimensional space, for every initial vector x(0) G X, the sequence {x(?t,)} 
generated by x(n + 1) = T(x(n)) converges to x* . 

As can be seen from Proposition |2l para-contractivity does not yield any estimate of the rate of 
convergence to the fixed point. 

B. Standard interference functions 

The standard interference functions were introduced by Yates ["4] to study various extensions 
of the basic distributed power control problem. 

Definition 1 (Standard Interference Function f4l) A function I : — )■ is called a 
standard interference function, if for all p > the following properties are satisfied: 

• Positivity: I(p) > . 

• Monotonicity : If p > p', then I(p) > I(p') . 

• Scalability: For all a > 1, al(p) > I(a!p). 
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The main convergence result for standard interference functions can be summarized as follows: 

Proposition 3 ( [4]|) Let I be a standard interference function and consider the iteration 

p(n + l) = l(p(n)). (2) 

Then, if (O has a fixed-point, this fixed-point is unique and the iterates {p(n)} produced by ^ 
converge to the fixed-point from any initial vector p(0). 

Note that contrary to the result for contraction mappings, the existence of fixed-points has to 
be verified separately, and no guarantees about the convergence rate of the iterates to the fixed- 
point are given. Already this should raise the suspicion that standard interference functions do 
not define contraction mappings. The following simple example establishes that this suspicion 
is indeed correct. 

Example 1 (Standard interference functions are not contractive) Consider I(p) = 2p + 1 
where 1 is the vector with all components equal to 1. This is a standard interference function, 
but ||I(p) — I(p')ll = 2||p — p'||, so it is neither contractive nor para-contractive. 

The interference function in Example [H does not have a fixed-point in the positive orthant, 
hence does not contradict Proposition [3l The next example shows that even a standard interference 
function has a positive fixed point, the iteration ^ may converge at a sub-linear rate. 

Example 2 Consider I{p) = -^^^^\_2) [0, oo). I(p) is a standard interference function since 
it satisfies positivity, scalability, and monotonicity. I {2) = 2, so p* = 2 is a fixed point. It is easy 
to verify that I'{2) = 1, I" (2) = 0, and I^^^2) = Therefore, the iteration Q converges 
sub-linearly to p* / fiJl Theorem 4]. 

Example [21 provides motivations for seeking stronger conditions than standard interference func- 
tion to ensure contractivity, hence linear convergence rates of the iterations. To this end, one 
could certainly make a separate analysis of contractivity of the particular interference functions at 
hand. However, if one can prove contractivity, particularly in the weighted maximum norm, then 
the interference function framework brings little additional value. The beauty of the framework 
lies in the easily verifiable conditions that guarantee synchronous and asynchronous convergence. 
Next, we will show that a slight reformulation of the scalability condition ensures contractivity. 
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III. Contractive interference functions 
We propose to study a class of interference functions which we call contractive. 

Definition 2 A function I : M.^ — )• M.^ is said to be a contractive interference function if it, for 
all p > 0, satisfies the following conditions: 

• Positivity: I(p) > . 

• Monotonicity: If p > p', then I(p) > I(p') . 

• Contractivity: There exists a constant c G [0,1) and a vector v > such that for all e > 0, 
I(p + ev) < I(p) + cev. 

Note that the two first conditions are the same as for standard interference functions, but the 
scalability condition has now been replaced by contractivity. The following example shows that 
contractivity, in general, does not imply scalability. 

Example 3 Consider the function 



Hp) 



5P ifi "f" inn 1 A — 



2^ 16 '100' 4 

which is contractive (c = |). However, the scalability property does not hold for a = 2 and 
p = I since I{ap) <^ al{p). 

However, with the additional requirement that : M.^ — t- IR+ is concave for all i = 1, . . . , K, 
the contractive interference function is also standard. This result follows immediately from the 
fact that positivity and concavity imply scalability [fT6ll . 

As shown in the next theorem, contractive interference functions define contraction mappings, 
which implies that the associated iterations ^ have unique fixed-points and linear convergence 
rates. 

Theorem 1 If I is a contractive interference function, then it has a unique fixed point p* and 
for every initial vector p(0), the sequence p{n + 1) = l(p(n)) converges linearly to p* as 
follows \\p{n) - pir^ < c"||p(0) - pII^. 

Proof: Let p ^ p' . Since \pi — p'i\ < ||p — p'H^fi for all i, we have 

P = P+P-P' 

^ ' I II 'liv 

< p + Hp - p 11^ V . 
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Note that ||p — p'||^ > 0. The mono tonicity and contractivity properties imply 

i(p) < i(p'+i|p-p'ir^v) 
< i(po+c||p-p'iro,v. 

By interchanging the roles of p and p', I(p') < I(p) + c||p — p'||^ v. So for all components of 
I(p),wehave \Ii{p) - Ii{p')\ < c||p-p'||^v, which implies that ||I(p')-I(p)||L < cUp-p'H^o- 
Therefore, according to Proposition [H / is contractive and the result follows. ■ 

To emphasize the modulus c of the contraction mapping, we say that a function is a c- 
contractive interference function. Note that convergence rate is directly related to the number of 
iterations required for the algorithm to converge. Specifically, if we define the convergence time 
of the iteration as the the smallest t such that ||p(t)— p*||^ < 5, then, with Rq = ||p(0)— p*||^, 
||p(t) — p*||^ < c*i?o- So if c < 1, < In We can see that the convergence time goes 
to oo as c tends to one. 

To show that the concept of contractive interference functions is useful, we will now show 
that it is readily applied to several distributed power control algorithms from the literature. We 
assume K users and R base stations and a common radio channel. 

First, consider fixed assignment interference functions 

T f \ SjV* ^raPj +Vri . 

where is the user z's base station, Grj is the link gain between base r and user j, is the 
target SINR of user i, and 77,. is the background noise at base r. Under fixed assignment, each 
user i has an assigned base station G {1, . . . , R}. Equation ([3]) can be rewritten as 

K 



where N^'''^ = -fi^, and 



j 7^ h 



0, 3 = ^- 



(r) 

Define M as an K x K matrix that has M-j as its elements. We then have the following result. 
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Theorem 2 Tjf ||M||^ < 1 for some v > 0, then linear interference functions ^ are c- 
contractive interference functions with c= ||M||^. 

Proof: It is clear that linear interference functions are positive and monotone. Furthermore, 

i=i i=i 
< /,(p) + e||M||>, , 

where we used Proposition |4l^b) in the appendix to get the last inequality. Hence, they are also 
contractive with c = ||M||^. ■ 
Since matrix M is a nonnegative square matrix, the following Proposition holds. 



Proposition 4 ( hl3\ Ch. 2]) Let A eW^^'^ be a nonnegative matrix. Then 

(a) There exists a positive vector v such that \\A\\1^ < 1 if and only if p{A) < 1. 

(b) Let V be a positive vector Then X]j=i^«j^j — II^IlL^* /'^'" 

According to Proposition |4l p(M) < 1 is a necessary and sufficient condition for the existence 
of a positive vector v for which ||M||^ < 1. Any vector of the form (/ — M)^^x where x > 
can then be chosen to be v [17, Proposition 1]. If the matrix M is irreducible, which is often a 
reasonable assumption (since we are not considering totally isolated groups of links that do not 
interact with each other), it is worth noticing that the positive right Perron-Frobenius eigenvector 
V is such that p(M) = ||M||^ ( fT3^, Proposition 6.6]). In either case. Theorem [2] confirms that if 
p(M) < 1, then the fixed assignment iteration has a unique fixed point and a linear convergence 
rate. That is, there exists a set of powers to which all transmitters converge exponentially fast, 
such that all transmitters meet their QoS requirements (minimum SINR for successful reception). 
This also coincides with previous tailor-made analyses (for example d, 0, ifTSl ). 

The convergence rate that can be guaranteed using Theorem [T] depends on the choice of v, 
and hence on the norm in which we require the iterations to be contractive. In many cases, 
sufficient but more easily verifiable conditions can be derived by considering v = 1. For the 
linear interference function iterations, v = 1 yields the condition that ||M||oo < 1, i-C-, 

^ > % Vz = l,2,...K 

Since p(M) < ||M||oo, this condition is more conservative, but it is easily verifiable as we only 
need to check if each of the row sums of matrix M are less than 1. 
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We now verify that the minimum power assignment and macro-diversity interference functions 
from [j4| are also contractive under appropriate assumptions. 

• Minimum Power Assignment. At each step of this iterative algorithm, user i is assigned to 
the base station r at which its transmitted power is minimized. In this case we have 

/,(p)= min |x:<v,+ivr 

re{i,...,R} 

where M^J^ and N-^^ are defined as (H)). Define matrices Mi, . . . , Mr where the 

(r) 

element of Mr is equal to M-j . Mr is the normalized gain matrix when all users are 
assigned to base station r. Let v > and c = maXr{||Mr||^}. Then 

{K K 
j2Ml;^p^ + Nt^ + Y,Ml; 
i=i i=i 

< mm I J2 M-pPj + iVf ^ + e\\M, 

[ j=i 

< /i(p) + cevi . 

At this point we have shown that if a common positive vector v exists such that ||Mr||^ < 1 
for all r, then the interference function is contractive. The question now would be how to 
check for the existence of such a v > 0, and this will be addressed now. Let L be a set 
containing all possible mappings / : {1, . . . , K} — {!,...,/?}. Each mapping / represents 
the allocation of users to base stations and denotes that under mapping /, user i is 
assigned to base station We can construct the normalized gain matrix under assignment 
/, M', by letting the i*'^ row of equalt the z*^ row of M;(j). The following lemma provides 
a test for the existence of a common v > 0. 

Lemma 1 Let be the normalized gain matrix under assignment I E L. If p{M^) < 1 
for all I, then there exists a positive vector v such that ||Mr||^ < 1 for all r. 

Proof: Since is a nonnegative matrix with zero diagonal entries, then if p{M^) < 1 
for all /, all eigenvalues of are strictly inside the unit circle. Thus, the real part of 
Aj(M' — /) is negative for all i, and hence — / is Hurwitz. On the other hand, the 
off-diagonal elements of M^ — I are nonnegative; therefore, — / is also Metzler. As a 
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result, there exists a positive vector v such that (M^ — /)v < for all r |fT9l Theorem 4]. 
Therefore, M^v < v holds, which implies that ||M,. ||^ < 1 for all r. ■ 
By Lemma[Il contractivity of minimum power assignment iteration boils down to computing 
the spectral radius of matrices MK A more conservative condition can be derived by 
considering v = 1. Then, / is contractive if ||M,. ||oo < 1 for all r. In this case, one is only 
required to calculate the row sums of R matrices to establish contractivity. 
Macro-Diversity: The interference function is defined as 

1 



Z^r=l 



Similar to minimum power assignment iteration, we can define R matrices Mi . . . , Mr. Let 

V > and c = maXr{||Mr||^}. We have 



Z^7 = l ^J 



/,(p + ev)-/,(p) = 



< 



^"^^ (EjLiM<;>p,+ivf>)(EEiMij'fe+^".)+^r') 



Thus, if c < 1 for some v > 0, then the interference function is contractive. By Lemma [H 
p(M') < 1 for all / G L guarantees contractivity. 

While the convergence rate estimates for the interference functions that we have considered 
earlier are tight or appear to be tight, the situation for the macro diversity results is less clear. 
Assuming that a user's power contributes to its own interference, Hanly [20] showed that the 
macro-diversity interference function has a unique fixed point if J2f=i 7i < ^- This condition 
is insensitive to the channel gains and hence to the position of base stations. While we are 
unaware of any convergence rate estimates for the original or modified macro diversity 
interference functions, a related interference function has been investigated by Rodriguez 
et. al [I2TII . To understand the relationship with our results, let mr,i(p) = Yl,j^i^rjPj and 
re-write the macro-diversity interference function as 
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By introducing mi(p) = max,{mr^j(p)} and rj = max,!?],}, 

- mi(p) + ?7 

is a strict overestimate of the original macro diversity interference function. For this inter- 
ference function, some results have been obtained in ETI . It is readily verified that /i(p) 
is contractive, since for any e > 0, we have 

max, l^jVi GrjiPj + evj)^ - max, j^jyi ^rjPj^ 



/i(p + ev) - Ji(p) = 7i 



< 



£7: 



max, GrjVj^ 



This inequality can be rewritten as 

/.(p + ev) - 7,(p) < emax 1^ H^^^j \ , (5) 



where 



\ 0, j=z. 

Define matrix Hr with {i,jy^ element equal to H-j\ Let c = max,{||i7,||^}, then In- 
equality ([5]) becomes /j(p + ev) — Ii{p) < cevi. Hence, if p{H^^^) < 1 for all I G L, then 
by Lemma [U < 1 for all r and the interference function is contractive. The result 

derived in [|2T| is 7t ^ ' < 1 for all i,r, which is equivalent to ||//, ||oo < 1 for all r. 



This indicates that the stability condition in pTI is more conservative than our result. 
To show that our framework allows to go beyond the known results, consider the utility-based 
power control (UBPC) from [ilOil . The associated interference function is 



where ai is a price coefficient and fi^{x) is a decreasing function on [K_^,Ki\ for all i given 
by fi{SIRi) = U-{SIRi) where Ui is a utility function of user i. In their paper, Xiao et al. use 
a sigmoidal utility function 
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where bi = % — a^^ ln(aj7i — 1). Let M^J-* be defined as (H)) and let Mh be a matrix with {i^jY^ 

(r) 

element equal to biM-j'. We will next show that the framework of contractive interference 
functions will allow us to analytically bound the convergence rate, which has an immediate use 
for tuning the algorithm parameters. Specifically, we have the following result. 

Theorem 3 Consider the interference function /" defined in ©-(IT]). If c = p(Mb) < 1, then 
/" is a c-contractive interference function. 

Proof: Let v > 0. For all e > we have 



/ K \ ( / K 



E<p, + A^r)j 



= /r (p) + 6 \^ M^v,]^ (^^ M^v, + A^f ' j j , (8) 

where the first inequality comes from the fact that f,^^{x) is a decreasing function. Under 
the sigmoidal utility function ([7]), the maximum value of Ul{SIRi) occurs at point bi, so 
max/j~^(x) = argmaxa; = h. Thus, Inequality dS]) becomes 



/r(p + ev)</r(p)+e|^fj6,M,,^,j 



</r(p)+e||M,r^. 

Hence, if llAf^H^ < 1, then UBPC iteration linearly converges to a unique fixed point and the 
result follows. ■ 
The following numerical example illustrates our result. 



Example 4 We consider UBPC under the simulation scenario described in /TO'/. Here, four 
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5 10 15 20 

Number of iterations 



Fig. 1. Comparison of upper bound on convergence rate of UBPC obtained in Theorem [5] and the actual convergence rate of 
UBPC for the scenario considered in example 4. 



mobiles share a channel with link gain matrix G given by 



10-* 6.82 X 10-'^ 3.57 x lO^^ 2.12 x 10"^ 

1.52 X 10-^ 6.25 X 10-'' 3.51 x IQ-^ i.gg x 10"^ 

G = 

7.67x10-9 2.44x10-** 1.23x10^6 5.16x10-9 

2.63 X 10-^ 4.82 x 10-** 2.56 x IQ-'^ 3.28 x IQ-^ 

The noise power is 0.5 and the target SIRs of the users are 6, 6, 8 and 10 dB, respectively. We 
assume that four users in the system use sigmoidal utility function with parameters 1.02, 1.32, 
0.88 and 1.05, respectively. The price coefficient a is equal to 5000 /or all users. Fig. \l\shows 
a comparison of the norm between power vector and optimal power vector versus number of 
iterations for UBPC algorithm and the theoretical bound obtained from Theorem\3\ It can easily 
be seen that the distance between the iterates of UBPC algorithm and optimal power allocation 
at iteration n, ||p(n) — p*||^, is always less than c"||p(0) — p*||^ and decays exponentially. 



Another useful result shows that imposing an upper and lower bound on a contractive inter- 
ference function does not change the contractivity properties. 

Theorem 4 If I is a c-contractive interference function, then so is 



F(p) = max{pmin,min{pmax,I(p)}}- 
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Proof: I''(p) satisfies positivity and monotonicity. It remains to show that I^(p) satisfies 
contractivity. Let e > 0. The contractivity condition of I(p) implies 

F(p + ev) = max{pmin,min{pmax,I(p + ev)}} 

< max{pmin,min{p 

max) 

I(p) + cev}} 

< max{pmin, min{pmax, I(p)} + cev} 

< maxjpinin, min{pmax, I(p)}} + cev . 

Therefore, I9(p + ev) < I9(p) + cev. ■ 
Finally, we consider the distributed robust power control algorithm (DRPC) from [|22l . Denote 
the normalized channel gain between user i's base station and user j as M^^^ = M-^^ + AM-j\ 
where M-^^ is the nominal value, and AM^p is the perturbation associated with M -^^ The 
uncertainty set of Ai is 

M = {M\MieMi, i = l,...,K}, 

where Mi is the i^^ row of matrix M = [M^p], and Aii is the uncertainty set of M^. The 
interference function of DRPC is defined by 

Ii{p)= sup (Mfp)+Ar,, i = l,...,K. 

The convergence rate of DRPC algorithm is verified by the following theorem. 

Theorem 5 If c = sup^^^^ II ^11^ < 1 far some v > 0, then DRPC interference function is 
c-contractive. 

Proof: Similar to Theorem [21 ■ 
Note that if Ai is bounded and sup^^g^^^ p(M) < 1, then there exists a v > for which 
supjv^g^ Mv < V Il22l Theorem 1]. Thus, Theorem [5] coincides with the stability condition 
derived in (221. 

IV. Interference functions and para-contractions 

We have already shown that standard interference functions do not define contraction map- 
pings. However, the convergence results for standard interference functions are identical to those 
of para-contractions in Proposition |2l so there should be a link between the two. This section 
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shows one such link. In particular, we will demonstrate that a logarithmic change of variables 
s = ln(p) makes the iterations para-contracting in the new variables. Interestingly, this is the 
same change-of-variables that has been very useful in convexifying various resource allocation 
problems for interference-limited wireless systems [|TT]| . [fT2l . 

Theorem 6 Suppose that the interference function I is standard. Then the change of variables 
s = ln(p) and I(s) = ln(l(e^)) transforms the interference function I : p ^ I(p) into a 
para-contracting function / : s — )■ I(s). 

Proof: We rewrite the properties of standard interference function in the new coordinates. 
Al If Si < S2, then e^^ < e^'\ By monotonicity property of I(p) 



I(e^^) < I(e^2) ln(l(e^i)) < ln(l(e^2)) 



^ I(si)<I(s2). 



A2 For any e > and by using scalability property of I(p), we have 



I(s + el) 



< 



ln(l(e^e^)) 
ln(e^I (e")) 



el + I(s). 



For any Si, S2 G M", we have 



Sl < S2 + 



S2||ool- 



Let Sl 7^ S2. By properties A1-A2 of I(s) 



I(Si) < I(S2 



+ — S2II00I) 



< 1(82) + ||Si - S2||ool- 



By interchanging the roles of Si and S2 in the preceding inequality, we obtain 



1(82) < 1(81) + ||8i - S2||ool- 



Consequently for all components of I(s), we can write 



/i(8i) - /i(82)| < ||Si - 82 



00 ) 



I = 
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which implies that ||I(si) — I(s2)||oo < — S2||oo- Therefore, I(s) is para-contracting. ■ 
Theorem [6] helps us to understand that interference functions are para-contractions with respect 
to a certain metric. Specifically, we note the following: 

Corollary 1 Standard interference functions are para-contractions with respect to the metric 

(ic(p, p') = max 

i 

It is important to note that standard interference functions are para-contractions on the metric 
space induced by dc, irrespective of whether they have a fixed-point or not. To guarantee 
convergence of the iterates, we must verify that the iteration has fixed-points. The following 
theorem can then be useful. 

Theorem 7 Given a standard interference function I, if there exists a p' such that I(p') < p', 
then a fixed point exists. 

Proof: Let ^ = {p : < p < p'}. By the positivity and monotonicity properties, we 
have < I(p) < p'. Hence, / maps X to itself. Let us define function / : X i— )■ [0, oo] 
by /(p) = '^c(p,I(p))- Since / is continuous and X is compact, / attains its minimum. Let 
/(p) = minpgx /(p)- If P is not a fixed point of /, then by para-contractivity 

/(I(p)) = 4(l(p),l(l(p)) 

< 4(p,i(p)) = /(p). 

contradicting the fact that p minimizes /. Hence I(p) = p. ■ 

V. Two-sided interference functions 

Sung and Leung [[5l present a new class of functions called two-sided scalable interference 
functions which generalizes the standard interference functions to allow for simple and powerful 
analysis of certain opportunistic power control laws: 

Definition 3 (Two-sided scalable interference functions (H) AfunctionI : — > is called 
a two-sided scalable interference function, if for all p > 0, I(p) satisfies: 

• Positivity: I(p) > 0. 

• Two-sided scalability: For all a> 1, ^p < p' < ap implies ^I(p) < I(p') < al(p). 



In 
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Note how monotonicity and scalability conditions of standard interference function have been 
replaced by the two-sided scalability condition. However, every standard interference function 
is also two-sided scalable. The key convergence result reads as follows: 

Proposition 5 ( [51) Let I be a two-sided scalable interference function and consider the iter- 
ation If the iteration has a fixed-point p*, then this fixed-point is unique and the sequence 
{p(n)} generated by the iteration converges to for every initial value p(0). 

The convergence conditions for two-sided scalable interference functions coincide with those of 
para-contractions. In particular, neither existence of fixed-points nor any convergence rates are 
guaranteed. To guarantee existence and uniqueness of fixed-points along with convergence rates, 
the natural concept is to consider two-sided contractive interference functions. 

Definition 4 A function I : Mf^ — )■ is called two-sided contractive interference function if 
it, for all p > satisfies 

• Positivity: I(p) > 0. 

• Two-sided contractivity: There exists a constant c G [0, 1) and a vector v > such that for 
all e > 0, p' — ev < p < p' + ev implies that I(p') — cev < I(p) < I(p') + cev. 

The associated convergence theorem now reads 

Theorem 8 If I is a two-sided contractive interference function, then (O has a unique fixed 
point p* and the sequence {p(n)} generated by the iteration ^ converges linearly to p* from 
every initial value p(0). 

The proof follows similarly to the contractive interference function proof and is omitted in this 
paper. Next, we show how two-sided scalable functions relate to para-contractions. 

Corollary 2 Suppose that / : p — )■ I(p) is two-sided scalable. Then I is a para-contraction 
with respect to the metric dc{p, p') = maxj | In p-| . 

Our results in this section are related to the work by MoUer and Jonsson [|23l , who studied 
stability of higher-order power control laws. They demonstrated that in logarithmic variables, 
two-sided scalability implies global Lipschitz continuity of the interference function, and an 
alternative restriction allows to establish convergence rates and uniqueness of fixed-points. 
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VI. Asynchronous Power Control 



So far, we have examined synchronous power control algorithms. In this case, every component 
of the vector p is updated at every time step, using information of the transmit powers used 
by all transmitters in the previous iteration. However, a nice feature of the standard interference 
functions is that they also guarantee convergence in the absence of synchronization. In this 
section, we will demonstrate that contractive interference functions also converge asynchronously. 

Asynchronous computation models may be divided into totally asynchronous and partially 
asynchronous jTSl Chapter 6-7]. Let T be the set of times when some transmitter updates 
its power, and let T* C T be the times when transmitter i executes an update. To model that 
the transmitter might need to update its power using old information from other transmitters, 
let rj(t) be the time at which the most recent version of pj available to node i at time t was 
computed. Node i executes the update 



Definition 5 (Total Asynchronism fT3] Chapter 6]) The iteration Q z^' totally asynchronous if 
the sets T* are infinite for all i, and if {tk} is a sequence of elements ofT^ that tends to infinity, 



Loosely speaking, this assumption guarantees that no transmitter ceases to update its power, and 
that such updates eventually propagate to all other transmitters in the network. Yates showed 
that if an iteration involving standard interference functions converges synchronously, it also 
converges when it is executed totally asynchronously. A similar result holds for contracting 
interference functions: 

Theorem 9 Let I be a contractive interference function. Then, the iterates produced by ^ 
converge to the unique fixed-point under total asynchronism. 

This result is proven by noticing that contractive interference functions define max-norm 
contractions which converge under total asynchronism {e.g., [[T3l Page 434]). A similar result 
can be established for two-sided contractive interference functions using the same arguments. 

While the totally asynchronous convergence result is comforting, it does not quantify the longer 
convergence times one could expect with increasing information delays. An additional advantage 




(9) 



then it also holds that lim^^oo '^j{tk) = oo for all j. 
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of the contractive interference function framework is that such results can be developed. Note that 
in the totally asynchronous model, the delays can become unbounded as t increases. Consider, in 
contrast, a situation where all mobiles update their powers at each iteration but the information 
delay is guaranteed to be bounded. In particular, there exists a positive integer D such that 
t — D < Tj(t) < t for all i and j. The following theorem gives a bound on the convergence rate 
of contractive interference functions under this model of asynchronicity. 

Theorem 10 If I : M.^ is c-contractive, then from any initial power vector p(0), the 

asynchronous power control algorithm satisfies: 



ipH-Piroo<ciip(o)-p'^"^ 



1 



where c = c^+i. 



Proof: Contractivity of / implies that ^^^^\ < cmaxjgji ... \ElJ!A^ i = 1, . . . , K. Let 
di be a maximum communication delay between user i and other users, then the convergence 
rate of the asynchronous algorithm Q is the unique solution of the equation [13, pp. 442] 

p = max{cp''^\- ■ ■ , cp""'^'}. (10) 

Since D = maxjdi, ■ ■ ■ , dx}. Equation (fTOl) can be written as p = cp^^ and by letting c = p, 
the proof is complete. ■ 
Note that contractive interference functions converge under arbitrary bounded delays, and that our 
result provides an explicit bound on the impact that an increasing delay has on the convergence 
rate. Moreover, computing the convergence time for the asynchronous update, we find that 

Ts<{D+l)-^\n^. 

Hence, a communication delay of D results in a convergence time that is no more than D + 1 
longer than that of the ideal synchronous iteration. 

Wu et al [[22l Theorem 2], using the definition for the spectral radius, derive the convergence 
rate of the robust interference function they consider, when updated every K steps (corresponding 
to D + 1 steps in our setup). More specifically, the convergence rate E for the case with no 
uncertainties, is given by 
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which is equivalent to our derived convergence rate. When, however, there exists reduced message 
passing and increased robustness consideration, this is incorporated in matrix M, thus reducing 
the rate of convergence. Therefore, their derived bounds on the convergence rates for the cases 
they consider, justify our analysis and constitute a worked example of our proposed framework. 

VII. Conclusions 

This paper has explored the connections between the standard interference function framework 
and the theory for fixed-point iterations. We have shown that standard interference functions 
do not define contraction mappings and introduced contractive interference functions, a slight 
modification of the standard interference functions, that guarantee existence and uniqueness of 
fixed-point along with linear convergence of iterates. We have demonstrated that several important 
distributed power control algorithms proposed in the literature are contractive and derived the 
associated convergence rates. In some cases, such as linear iterations, the convergence rate 
coincides with known results from the literature that has been obtained using a detailed and 
tailored analysis. In other cases, such as the utility-based power control, we provide the first 
convergence rate estimates in the literature. This paper also provided a link between standard 
interference functions and para-contractions. This link involves a logarithmic change of variables 
(alternatively, analysing the iterations with respect to a specific metric), which coincides with 
the change-of-variables that has been so successful in convexifying resource allocation problems 
in interference-limited wireless systems. Associated results for two-sided scalable interference 
functions have also been given. 
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